NPS ARCHIVE 
1997 . 0 ^ 
AKIN, I. 



NAVAL POSTGRADUATE SCHOOL 
Monterey, California 




THESIS 



A ROBUST SYMMETRICAL NUMBER SYSTEM 
WITH GRAY CODE PROPERTIES 
FOR APPLICATIONS IN SIGNAL PROCESSING 



Thesis 

A3336 



by 

Ilker Aydin Akin 
September, 1997 



Thesis Advisor: 



Phillip E. Pace 



Approved for public release; distribution is unlimited. 




SwA5S^|ESCH00t 






REPORT DOCUMENTATION PAGE 



Form Approved OMB No 0704-0188 



Public reporting burden for this collection of information is estimated to average 1 hour per response, including the ume tor reviewing instruction, searching existing data sources, 
gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of this 
collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for Information Operations and Repons, 1 215 
Jefferson Davis Highway, Suite 1204, Arlington, VA 22202-4302, and to the Office of Management and Budget, Paperwork Reduction Project (0704-0188) Washington DC 
20503. 



1. AGENCY USE ONLY (Leave blank) 2. REPORT DATE 

September 1997 


3. REPORT TYPE AN D DATES COVERED 
Master’s Thesis 


4. TITLE AND SUBTITLE 

A ROBUST SYMMETRICAL NUMBER SYSTEM WITH GRAY CODE 
PROPERTIES FOR APPLICATIONS IN SIGNAL PROCESSING 


5. FUNDING NUMBERS 


6. AUTHOR(S) Akin, Ilker A. 


7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 
Naval Postgraduate School 
Monterey, CA 93943-5000 


8. PERFORMING 

ORGANIZATION 
REPORT NUMBER 


9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 
Research and Development Office 


10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 



11. SUPPLEMENTARY NOTES The views expressed in this thesis are those of the author and do not reflect the official 



policy or position of the Department of Defense or the U.S. Government. 



12a. DISTRIBUTION/AVAILABILITY STATEMENT 


12b. 


DISTRIBUTION CODE 


Approved for public release; distribution is unlimited. 







1 3. ABSTRACT (maximum 200 words) 



A new symmetrical number system with applications in parallel signal processing is investigated. The “Robust 
Symmetrical Number System” (RSNS) is a modular system in which the integer values within each modulus, when 
considered together, change one at a time at the next position (Gray code properties). Although the observed dynamic range 
of the RSNS is somewhat less than the optimum symmetrical number system, the Gray code properties make it particularly 
attractive for folding analog-to-digital converters. With the RSNS, the encoding errors (due to comparator thresholds not 
being crossed simultaneously) are eliminated, as is the need for the corresponding interpolation signal processing (reduced 
complexity). Computer generated data is used to help determine the properties of the RSNS. These properties include the 
largest dynamic range (number of distinct consecutive vectors), and the position of the largest dynamic range within the 
system. The position of the maximum unambiguous dynamic range is also quantified. Least squares analysis of 2 and 3 
moduli systems are used to formulate closed-form expressions for the dynamic range. To compare the advantages of the 
RSNS with previously published results, the transfer function of a 3-channel RSNS folding analog-to-digital converter 
architecture (m 7 =3, m 2 = 4 , and m 3 = 5) is numerically evaluated using SPICE. 





14. SUBJECT TERMS Residue Number System, Robust Symmetrical Number System, 
Relatively Prime Numbers, Pairwise Relatively Prime Numbers, Moduli, Analog-to- 
Digital Converter, Dynamic Range, Point Dynamic Range, Fundamental Period, Point 
Index, Gray Code, Least Significant Bit 


15. NUMBER OF 
PAGES 
102 


16. PRICE CODE 


1 7. SECURITY CLASSIFICA- 
TION OF REPORT 
Unclassified 


18. SECURITY CLASSIFI- 
CATION OF THIS PAGE 
Unclassified 


1 9. SECURITY CLASSIFIC A- 
TION OF ABSTRACT 
Unclassified 


20. LIMITATION OF 
ABSTRACT 
UL 



NSN 7540-01-280-5500 Standard Form 298 (Rev. 2-89) 



Prescribed by ANSI Std. 239-18 298-102 



I 



ii 



Approved for public release; distribution is unlimited. 



A ROBUST SYMMETRICAL NUMBER SYSTEM 
WITH GRAY CODE PROPERTIES 
FOR APPLICATIONS IN SIGNAL PROCESSING 



Ilker Aydin Akin 

Lieutenant Junior Grade, Turkish Navy 
B.S., Turkish Naval Academy, 1991 



Submitted in partial fulfillment of the 
requirements for the degree of 



MASTER OF SCIENCE IN SYSTEMS ENGINEERING 



from the 



NAVAL POSTGRADUATE SCHOOL 
September 1997 



JUL/lcY *NOX LIBRARY 

^cd2S TGRADUATE school 

^OK'Tcpcv CA t W943 ^ 101 

ABSTRACT 

A new symmetrical number system with applications in parallel signal processing is 
investigated. The “Robust Symmetrical Number System” (RSNS) is a modular system in 
which the integer values within each modulus, when considered together, change one at a 
time at the next position (Gray code properties). Although the observed dynamic range of 
the RSNS is somewhat less than the optimum symmetrical number system, the Gray code 
properties make it particularly attractive for folding analog-to-digital converters. With the 
RSNS, the encoding errors (due to comparator thresholds not being crossed 
simultaneously) are eliminated, as is the need for the corresponding interpolation signal 
processing (reduced complexity). Computer generated data is used to help determine the 
properties of the RSNS. These properties include the largest dynamic range (number of 
distinct consecutive vectors), and the position of the largest dynamic range within the 
system. The position of the maximum unambiguous dynamic range is also quantified. 

Least squares analysis of 2 and 3 moduli systems is used to formulate closed-form 
expressions for the dynamic range. To compare the advantages of the RSNS with 
previously published results, the transfer function of a 3-channel RSNS folding analog-to- 
digital converter architecture {mi =3, m 2 = 4, and m 3 = 5) is numerically evaluated using 
SPICE. 
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I. INTRODUCTION 



A. MODULAR NUMBER SYSTEMS 

This study addresses a new modular number system. Before going into the details 
of the new system, it is useful to explain the properties of the some of the other modular 
number systems and examine their applications in Signal Processing. One of the 
drawbacks of the modular signal processing is the introduction of large errors (encoding 
errors) in the decoded binary or decimal values when only small errors occur within each 
modulus. The encoding errors and their elimination are the main reason for defining a 
new modular number system. 

1. Residue Number System 

The residue number system (RNS) is representative of a sawtooth folding 
waveform with folding period equal to the modulus m. An introduction to the RNS has 
been given by Szabo [1]. It has been established as an important tool in parallel 
processing applications and high-speed computations. The RNS can serve as a source for 
parallel operations (moduli) that are of smaller computational complexity. Each parallel 
operation has a different modulus and requires only a precision in accordance with that 
modulus. The folding waveform of the RNS is generated as 

a = r (mod m) (1) 

where 0 < r < m, where m is the modulus (channel value) which are PRP and r is the 
residue of a modulo m. Figure 1.1 shows the RNS folding waveforms for mi = 4 and m 2 
= 5. Going from integer values within each modulus to decimal value requires the 
Chinese Remainder Theorem (CRT) [2-5]. The integers r within each channel (modulus) 
are representative of a folded waveform with the period 

P rns = m (2) 

The dynamic range (region of no ambiguities) of the RNS is 

N 

^ RNS ~ 1 1 m i 
i=\ 

th 

and m, is the i pairwise relatively prime (PRP) modulus (channel value), 1 <i <N. 



( 3 ) 



Sawtooth Folding Waveform Sawtooth Folding Waveform 



m j = 4 




m 2 = 5 




Figure 1.1: The Folding Waveform of the RNS for mi = 4 and m 2 = 5 
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2. Symmetrical Number System 

The SNS is also composed of a number of pairwise relatively prime (PRP) 
moduli. The integers within each SNS modulus however are derived from a 
symmetrically folded waveform [6]. The integer values within each SNS modulus are 
derived from a symmetrical folding waveform and therefore incongruent modulus m. 
Due to the presence of ambiguities, the set of integers within each SNS modulus does not 
form a complete residue system by themselves. It is well known that the inclusion of 
additional redundant moduli can effectively detect and correct errors that may occur 
within a RNS representation of a number. The SNS formulation is based on a similar 
concept, which allow the ambiguities that arise within the SNS to be resolved by using 
various arrangements of the SNS moduli. The integer values within each SNS folding 
waveform is given as 



x m % odd 



0 , 1 ,..., 



m 

2 



m 

2 



,..., 2,1 



, when m is odd 



m m 

0 , 1 ,...,-,-- 1 ,... , 2,1 



, when m is even 



(4) 

(5) 



where |_*J indicates the greatest integer less than or equal to x Both have size 1 x m 
and consist of the symmetrical remainder elements Xh, 0 < h < in. Figure 1.2 shows the 
SNS folding waveform for mi =4 and m 2 - 5. 

The integers within each modulus are representative of a symmetrically folded 
waveform with the period of the waveform equal to the modulus 

Psns = m ( 6 ) 



where m is the modulus (channel). 

The dynamic range of the SNS M SNS < Mosns is given as 
(i) For PRP moduli m„ with one of the moduli even, then 



M sns = min- 



m, 



9 Fk + Ik 

L 1=2 l=j+\ 



(7) 



where j ranges from / to N -J and nia, 111 / 3 , ... m,^ range over all permutations (2, 3 

NJ. The product of j, from 2 to l is empty and its value is I . 
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SNS Folding Waveform _ SNS Foldi ng Waveform 



( 8 ) 



(ii ) If all the PRP moduli in the system are odd, then 




where j ranges from 1 to N - I and mu, m mu, ... mm range over all permutations ( 1 , 2, 

3 NJ. 



m 



l 





m 2 = 5 




Figure 1.2: The Folding Waveform of the SNS for m t = 4 and m 2 = 5 
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3. Optimum Symmetrical Number System 

The optimum symmetrical number system (OSNS) formulation is a direct 
consequence of the need to extract the maximum amount of information from a 
symmetrically folded waveform [7], The OSNS scheme is also composed of a number of 
PRP moduli m. Integers within each OSNS modulus is defined as 

x m = [ 0, 1 m - 1, m- 1 , .... 1, 0]. (9) 

Figure 1.3 shows the OSNS folding waveform for mi = 4 and m 2 = 5. The integers 
within each SNS modulus are representative of a symmetrically folded waveform with the 
period of the waveform equal to the twice the modulus and 

PosNS = 2m (10) 

where m is the modulus. Due to the presence of ambiguities, the integers within (9) do 
not form a complete system of length 2m by themselves. The ambiguities that arise 
within the modulus are resolved by considering the paired values from all moduli 
together. By recombining the N moduli, the OSNS is rendered a complete system having 
a one-to-one correspondence with the residue number system. For N equal to the number 
of PRP moduli, the dynamic range of the OSNS is 

N 

M OSNS = n-, (11) 

1=1 

where m, is the i' th modulus, 1 < i <N and N is the number of moduli. The dynamic range 
is also the position of the first repetitive moduli vector. 

4. Application for Modular Number Systems in Signal Processing 

A modular number system such as the RNS has been used in high performance 
digital signal processing [8]. Applications include a high rate arithmetic processor, 
digital filters, memory compression, and image processing. Recently, symmetrical 
number systems have been investigated in such application as Analog-to-Digital 
Converters (ADCs) [7, 9-13], Digital Antennas [14], Undersampling Receivers [6, 15], 
and Direction Finding Antennas [14]. The symmetrical number systems offer a 
considerable amount of flexibility in implementation. For higher resolution, faster 
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OSNS Folding Waveform OSNS Folding Waveform 



conversion speed, and lower power dissipation the symmetrical number systems play an 
important role in folding ADCs applications. 



m j = 4 




m 2 = 5 




Figure 1.3: The Folding Waveform of the OSNS for my -4 and m? = 5 
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5. Encoding Errors of Modular Number Systems 

Each of the above modular number systems for which the integer values within 
each modulus, when considered together, change one or more at the next position is 
shown in Table 1.1. Consequently, when some integers at a position to change, do 
change, while others do not, there will be a large error (encoding error). Table 1.1 shows 
the integer values within the modular systems for moduli m t = 4 and m 2 = 5. 

When the input is 1, the integer values within the modular number systems are 1 
for both moduli. By increasing the input from 1 to 2, all integer values within the 
modular number systems change from 1 to 2. If one of the integer values does not 
change, the new moduli pairs, for example (1,2) or (2,1) correspond to new input values. 
For the RNS, the input values can be 6 or 17. For the SNS, the input values can be 6, 13, 
14, or 17. And for the OSNS, the input values can be 17 or 18. These values cause large 
errors (encoding errors) that must be dealt with. These large encoding errors present a 
number of difficulties when used for signal processing. The goal of this thesis is to 
develop a new number SNS formulation that minimizes the encoding errors. To examine 
the performance of the new number system formulation a folding ADC architecture is 
investigated. 

B. PRINCIPAL CONTRIBUTION 

A new symmetrical number system with applications in parallel signal processing 
is investigated. The “Robust Symmetrical Number System” (RSNS) is a modular system 
in which the integer values within each modulus, when considered together, change one 
at a time at the next position (Gray code properties). Although the observed dynamic 
range of the RSNS is somewhat less than the optimum symmetrical number system, the 
Gray code properties make it particularly attractive for folding analog-to-digital 
converters. With the RSNS, the encoding errors (due to comparator thresholds not being 
crossed simultaneously) are eliminated, as is the need for the corresponding interpolation 
signal processing (reduced complexity). Computer generated data is used to help 
determine the properties of the RSNS. These properties include the largest dynamic 
range (number of distinct consecutive vectors), and the position of the largest dynamic 
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range within the system. The position of the maximum unambiguous dynamic range is 
also quantified. Least squares analysis of 2 and 3 moduli systems is used to formulate 
closed-form expressions for the dynamic range. To compare the advantages of the RSNS 
with previously published results, the transfer function of a 3-channel RSNS folding 
analog-to-digital converter architecture (m/ =3, m 2 = 4, and m 3 = 5) is numerically 
evaluated using SPICE. 

C. THESIS OUTLINE 

This thesis is organized in six parts. Chapter II provides the definition of the 
robust symmetrical number system (RSNS). An expression for the fundamental period is 
given and proven. In this chapter, the algorithm, which finds the dynamic range, point 
dynamic range, and fundamental period is discussed. The beginning and ending of the 
discrete states without redundancy for any given number of channel and shift values of 
the RSNS, is presented. 

Chapter HI investigates the 2-channel (moduli) RSNS. A specific number 
combination is considered. The RSNS performance using different moduli, the same 
distance apart is evaluated to find a closed-form expression dynamic range in the general 
case. 

Chapter IV deals with the 3-channel (moduli) RSNS. A specific number 
combination is also considered in order to find a closed-form expression for the dynamic 
range. 

Chapter V is based on the simulation of a 3-channel RSNS ADC in SPICE. The 
simulation channel and the shift values are m/ = 3, sj = 0, m 2 = 4, S 2 = 1, and m 3 = 5, s 3 
= 2, where m shows the channel values and s shows the shift. The folding circuit and 
comparator stages are designed to fold and quantize the incoming signal in a RSNS 
format. 

Chapter VI summarizes the results and conclusions of this study and provides 
several recommendations for follow-up research. 
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Table 1.1: Modular Number Systems to Decimal Mapping for mi = 4 and m 2 = 5 
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II. ROBUST SYMMETRICAL NUMBER SYSTEMS 



A. DEFINING THE RSNS WAVEFORM 

In the robust symmetrical number system (RSNS), N different periodic 
symmetrical waveforms are used with, pairwise relatively prime (PRP) integer, lengths 
mi, m 2 , .... mu. The RSNS is based on the following sequence 

Xml = [ 0, 1,2 m - 1, m, m - 2, 1 ] (12) 

where jc m/ is the row vector and m is a positive integer (m > 0). In an V-channel RSNS, 

where N >2, the basic sequence for i' h channel, with modulus m is 

x m 2=[0, 0,..., 0, 0, 1, 1 1, 1,..., m, m m, m .... 1, 1 1, 1]. (13) 

In this sequence each value in the ;c m / row vector is put N times in succession. This 

sequence is repeated in both directions, forming a periodic sequence with the period 

P rsns = 2mN. ( 1 4 ) 

An W-channel RSNS is formed of vectors by picking N moduli (channels) m„ and N shift 
values Sj, 1 <i <N, where {si, S2,..., suf forms a complete set of residues modulo N. The 
V h vector 



a 



I h 



V h = 



(15) 



_ a m i_ 

th 

where is the (h - si) term in the periodic sequence for /»„ 1 <i <N. Since the i row is 
periodic with period 2Nm\, 1 <i < N, it is clear that if the length that is simultaneously a 
multiple of 2Nm /, 2Nm 2 ,..., and 2Nmu is taken, then the W-channel RSNS vector will 
repeat itself. Since the fundamental period for channel i is 2Nmi, 1 < i <N, it follows that 
a period for the RSNS vectors must be a multiple of each of 2Nm,. Therefore the 
fundamental period for the RSNS is 

PFrsns = [2m iN, 2m 2 N, 2m N N] ( 1 6) 
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where [nj, n 2 , ..., n n] is the least common multiple of «/, n 2 , .... n n. From number theory 
[2-5], 

PFrsns = 2N[ mi, m 2 , .... m N ]. ( 17 ) 



Considering a single channel the discrete states of the robust symmetrical number 
system are given as 



g(ni,s : ,n) 



r x — s i 

L N 

< 

2Mn. + N - n + S; - 1 



for \ + s t <n < Nni + .v 
for Nm t + 1 + s i < n < 2 Nm i + s i 



(18) 



where m, is the channel modulus and Sj is the corresponding shift for 

si = 0, 1, 2,..., N-l (modN) (19) 

and N >2 is the number of channels in the system. The values {si, s 2 ,..., s^} must form a 
complete residue system modulo N. Here n represents the normalized input voltage. 
Because of the relative property of the shifts, one of the shift values will be set equal to 0. 

Upon inspection of any 2-channel and 3-channel RSNS, it is seen that if V h (the 
vector at position h ) equals Vh +k , then k is a multiple of 4 and 6, respectively. These are 
particular examples of a general rule: 

Rule 2.1: Suppose that Vh, Vh+k are the vectors at position h and h + k, respectively, in an 
N-channel RSNS. IfVh = Vh+ k , then k is the multiple of2N. 

The proof is given below, but first recall that in an 77-channel RSNS each term in a 
given channel occurs consecutively N times. Consider the N consecutive, equal terms 
form from left to right as lying in the first, second,. N 1 ' 1 position. Since 2N is a multiple 
of N, this rule shows that at any redundancy the redundant term lies in the same relative 
position as the original term. That is, if Vh = Vh+k, 





- 




- - 




a \h 




a \,h+k 


II 


• 


V = 

’ v h + k 






_ a Nh _ 




^ a N,h + k _ 
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th 

and a At h occurs in the i position in row A. Then, a A .h = a A.h+k and a A .h+k also occurs in 
the i ,h position in row A, 1 <A <N. 

The proof of the Rule 2.1 : 

In the fundamental building block of the RSNS, the sequence by (12), the value at 
j (0 < j < m) equals to the value at (j + k) when j + k = m + (m- j), i.e. when k = 2m - 2j, 
an even number. This fundamental building block has length 2m, and it is extended 
periodically. The values 0, and m repeat every 2m terms, and no more. Other values 
repeat twice as often, but always an even number of terms apart. When the sequence is 
expanded by the /V-fold repetition of each term, then all terms that repeat in the same 
relative position are a multiple of 2N apart. 

Suppose that ah = k > 0, and that a* is in the first (left) position. Then «/,+* 
will be just i terms to the right of the term in the first position, where 0 < i < N. By the 
conclusion of the last paragraph k = 2Nj + i for some integers i, 0 <i < N, and j >0 ( and 
not both 0). Similarly, if a* = ah+k, where a* is in the N 1 ' 1 (right) position, then it follows 
that k = 2Nj - i , for some i ,0 <i < N, and j > 0. 

Now return to the full /V-channel problem. Suppose V/, = Vh+k with k > 0. By the 
definition of the RSNS one of the terms in (15) is in the first position, say a Ai h, and one is 
in the N ,h position, say as > Since Vh+k equals V*, a At h - a A ,h+k and as,h = The first 

of these equations implies that k = 2Nj + i for some i,0 <i < N and j > 0, and the second 
equation implies that k = 2Nj - i , for some i , 0 < i < N, and j > 0. Putting these 
together: Since 2Nj + i and 2Nj - i are both equal to k, 2Nj + i = 2Nj - i . Thus i + i = 
2N(j - j). However, 0 <i + i < 2N, so / + / = 0 and j = j . Since i and i are nonnegative 
both must be zero. Therefore k is a multiple of 2N, as claimed. 

The discrete states of the RSNS are indexed starting from the vector, which has all 
zero values. The index value starts from 0. These index values are put in a row matrix h, 
which lies in { ..., -2, -1, 0, 1, 2, 3,... } and these values are called the point index. For 
2-channel RSNS, all point indexes determine 2 non-negative integer numbers, which are 
the discrete states of the RSNS. Similarly, /V-channel RSNS point indexes determine N 
discrete states of the RSNS. 



13 



Figure 2. 1 shows the 2-channel RSNS folding waveforms and the discrete states 
within the RSNS. The thresholds shown on the vertical axis represent the integer value 
within each RSNS modulus. The integer values occur 2 times in succession. In this 
example the channel and corresponding shift values ar e mj = 4, si = 0 and m 2 = 5, S 2 = 1 . 

Figure 2.2 shows 3-channel RSNS folding waveforms and the discrete states 
within the RSNS. The thresholds shown on the vertical axis represent the integer value 
within each RSNS modulus. The integer values occur 3 times in succession. In this 
example the channel and corresponding shift values are mi = 3, Sj = 0, m 2 = 4, S 2 = 1 
and m3 = 5, S3 - 2. 

To demonstrate the performance of the 2 and 3-channel RSNS in detail, the 
RSNS computer algorithm is executed and Tables 2.1 and 2.2 show the distinct vector 
positions. Each row is considered separately. The first and second column of a given 
row give the beginning and ending point indexes containing consecutive sequences of 
distinct vectors of N states. Furthermore, these positions can not be extended without 
introducing an ambiguity. The total length of the consecutive sequence of distinct 
vectors, starting at the beginning points (the point dynamic range) is shown in the third 
column. Note the maximal number of distinct vectors in Table 2.1 is 21 and in Table 2.2 
is 43. 
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mi = 4, si =0 




m 2 = 5, s 2 = 1 




m=5 1001 1223344554433221 1001 1223 3 4 4(5=7) 
Figure 2.1: The RSNS folding waveforms for m/ = 4, Si = 0 and m 2 = 5, s 2 = 1 
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nti = 3, si = 0 




m 2 = 4,s 2 = 1 




m 3 = 5, s 3 -2 




m—4 10001 1 12223334443332221 1 10001 1 ( s=l ) 
m=5 11000111222333444 5 5 54443332221 (s=2) 
Figure 2.2: The RSNS folding waveforms for mi = 3, si = 0, m 2 = 4, s 2 = 1 and 

m 3 = 5, s 3 = 2 
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8 


28 


21 


14 


29 


16 


15 


33 


19 


31 


42 


12 



Table 2.1: The Position of Distinct Vectors for the 2-Channel Case RSNS 

(mi = 4, si = 0, and m 2 = 5, s 2 = 1 ) 



Beginning Point 
Index 


Ending Point 
Index 


Total 

Length 


2 


29 


28 


7 


30 


24 


8 


41 


34 


13 


51 


39 


23 


52 


30 


24 


57 


34 


41 


65 


25 


61 


103 


43 



Table 2.2: The Position of Distinct Vectors for the 3-Channel Case RSNS 

(mi = 3, si = 0, m 2 = 4, s 2 = 7, and m 3 = 5, S 3 = 2) 
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B. THE DYNAMIC RANGE OF AN RSNS 



The dynamic range of the RSNS is the maximum number of distinct vectors 

without a redundancy. The dynamic range of the RSNS is shown with the symbol M RSNS . 
The point dynamic range of the RSNS is the maximum number of consecutive vectors 
without a redundancy starting at a given vector (point index). The point dynamic range of 

the RSNS is shown with the symbol M P RSNS . 

The performance of a 2-channel RSNS is shown in Table 2.3 for m 2 = 4, sj = 0 
and m 2 = 5, s 2 = 1 . The selection of the shifts and their permutations among the 2-moduli 
has no effect on the dynamic range. The point index (beginning-ending) of the set of 
vectors that gives the dynamic range, however, is different. 



Moduli 

(md 


Shift 

to) 


Dynamic 

Range 

( M RSNS ) 


Fundamental 

Period 

(Prsns) 


Beginning - Ending Point 
Index of 
Dynamic Range 


4 


0 


21 


80 


5-25 


5 


1 








4 


1 


21 


80 


14-34 


5 


0 









Table 2.3: The Performance with the Permutation of Shifts for the 2-Channel Case 

RSNS 



The performance of a 3-channel RSNS is shown in Table 2.4 for mj = 3, Sj = 0, 
m 2 = 4, s 2 = 1 and m 2 = 5, s 2 = 2. The selection of the shifts and their permutations 
among the 3-moduli also has no effect on the dynamic range. The point index 
(beginning-ending) of the set of vectors that determines the dynamic range is different. 
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Moduli 


Shift 


Dynamic 


Fundamental 


Beginning - Ending Point 


(>H|) 


(Si) 


Range 


Period 


Index of 






( Mrs NS ) 


(Prsns) 


Dynamic Range 


3 


0 








4 


1 


43 


360 


61-103 


5 


2 








3 


0 








4 


2 


43 


360 


142-184 


5 


1 








3 


1 








4 


0 


43 


360 


8-50 


5 


2 








3 


1 








4 


2 


43 


360 


134-176 


5 


0 








3 


2 








4 


0 


43 


360 


180-222 


5 


1 








3 


2 








4 


1 


43 


360 


81-123 


5 


0 









Table 2.4: The Performance with the Permutation of Shifts for the 3-Channel Case 

RSNS 
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C. A WEAKNESS OF THE RSNS WITH RESPECT TO THE INCREASE OF 
THE NUMBER OF CHANNELS 

The integers fall into two classes, even and odd parities. The integer values for any 
number of RSNS channels examined by a parity point of view. For the 2-channel ideal 
case, their parities distribute 2 2 = 4 different ways. These are even-even, even-odd, odd- 
even, and odd-odd. The possible parity sets for the 2-channel RSNS are shown in Table 
2.5 and these parity sets repeat themselves. In the 2-channel RSNS, the efficiency of the 
parity sets is same as the ideal one. 



Index 


mi 


m 2 


1 


0 


0 


2 


0 


1 


3 


1 


1 


4 


1 


2 


5 


2 


2 


6 


2 


3 


7 


3 


3 


8 


3 


4 


9 


4 


4 


10 


4 


5 


1 1 






12 






13 








Table 2.5: Demonstration of the parity combinations for N = 2 with sy = 0 and s 2 = 1 

For the 3-channel ideal case, their parities distribute 2 3 = 8 different ways. These 
are even-even-even, even-even-odd, even-odd-odd, even-odd-even, odd-even-even, odd- 
even-odd, odd-odd-even, and odd-odd-odd. There are only 6 possible parity sets for the 
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3-channel RSNS (restricted parity combinations) and they repeat themselves as shown in 
Table 2.6. The parity sets even-odd-even and odd-even-odd will never show up in this 3- 
channel RSNS. In the 3- channel RSNS, the efficiency of the parity sets is 25% less than 
the ideal one. 



Index 


m, 


m 2 


m 3 


1 


0 


0 


0 


2 


0 


0 


1 


3 


0 


1 


1 


4 


1 


1 


1 


5 


1 


1 


2 


6 


1 


2 


2 


7 


2 


2 


2 


8 


2 


2 


3 


9 


2 


3 


3 


10 


3 


3 


3 


11 


3 


3 


4 


12 


3 


4 


4 


13 








14 








15 










Table 2.6: Demonstration of restricted parity combinations for N = 3 with 

si = 0, S2 = 1 and S3 = 2 



In the N-channel RSNS each term in the vector V^+^by (20) differs by exactly one 
from corresponding term in V h by (15). Therefore, each term in Vh+N has opposite parity 
from the corresponding term in V/,. It follows quickly that each term in Vh+ 2 N has exactly 
the same parity as the corresponding term in Vh. From this result, it can be said that the 
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parities of the terms in any V /, will be exactly the same as the parities of the corresponding 
terms in one of V], V 2 , ..., Vn- This conclusion ends the proof that there will be just 2N 
different parity sets represented, out of a possible 2 N for A-channel RSNS. 

D. RSNS COMPUTER ALGORITHM 

In this study, one of the important goals is to be able to find a closed-form 
expression for the dynamic range and fundamental period of any number of RSNS 
channels. The fundamental period of any number of RSNS is given in Section A. The 
RSNS computer algorithm (rsns.m) is constructed in order to examine the performance 
of the RSNS and find the dynamic range for any N > 2 RSNS systems. 

In rsns.m, the channel and shift values, which are going to be used must be entered 
and the discrete states of the RSNS are calculated by using these values along with (18). 
The calculated discrete states of the RSNS are put into a matrix. This matrix is extended 
two times the fundamental period of the RSNS channels. 

The next step in the algorithm is to find the point index of the discrete states within 
the extended matrix. After finding the point index, the data is put into a new matrix. The 
point index of each discrete state and its redundancy is assigned as the beginning and 
ending point index, respectively. 

The values of the new matrix are sorted and rearranged such, that the beginning and 
ending point index are monotonically increasing. Then, the length at each point index is 
found (Ending point index - Beginning point index + 1). The maximum value of these 
lengths is the dynamic range of the RSNS. 

The algorithm rsns.m also displays the dynamic range, fundamental period, entered 
channel and shift values, the beginning-ending point index values of the discrete states 
without redundancy and the point dynamic range of the RSNS. The program structure 
and listing of the rsns.m is shown in Appendix. Tables 2.1 and 2.2 are the example runs 
of the rsns.m algorithm. 
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III. 2-CHANNEL ROBUST SYMMETRICAL NUMBER SYSTEM 



A. SPECIAL CASE: MODULI SPACING 1 OR 2 APART 

In this chapter, the 2-channel RSNS is addressed. The 2 moduli are RP numbers 
and the different pairs are chosen to have the same distance between each other. The 
algorithm rsns.m is used to find the dynamic range. When the moduli have a difference 
( A ) of 1 or 2, the following computer results have been observed as shown in Table 3.1 
and 3.2. 



m, 


m 2 


Dynamic Range 


1 st Difference 


3 


4 


15 


6 


4 


5 


21 


6 


5 


6 


27 


6 


6 


7 


33 


6 




Mrsns = 3m < 


+ 3m 2 - 6 





Table 3.1: The dynamic range of mj and m 2 = mi + 1 where m t > 3 (A = 1) 



m 2 


mi 


Dynamic Range 


1 st Difference 


7 


5 


29 


12 


9 


7 


41 


12 


11 


9 


53 


12 


13 


11 


65 


12 


Mrsns - 3 m i + 3m 2 - 7 



Table 3.2: The dynamic range of mi and m 2 = m t + 2 where mj >5 (A = 2) 



These relationships show a uniform increase of 3 with each unit increase of value 
in mj and m 2 . The formulas 
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= 3m i + 3 m 2 - 6 



( 21 ) 



Mr SNS 

Mrsns ~ 3mi + 3m 2 - 7 (22) 

will work for the first and second sets of cases which are shown in Table 3.1 and 3.2, 
respectively. These cases can be expressed in terms of the lower value, mj and a 
difference, A = m 2 - mi. In the first set A = 7, and in the second set A = 2. Thus, (21) and 
(22) can be written as 

Mrsns - 3 }n i + 3(mi + A) - 6 = 6m 1 - 3 (23) 

Mrsns ~ 3 m i + 3(mj + A) - 7 = 6mi - 1. (24) 

These can be combined as 

Mrsns = 6 m i + 2 A- 5 (25) 

and now A can be replaced with m 2 - m 2 , yielding 

Mrsns ~ 6 m i + 2(m 2 -m,) - 5 = 4m 1 + 2m 2 - 5 (26) 

and this formula works for both A = 1 and A = 2. This formula is checked against the 
computer generated results. For example, when mi=100 and m 2 =101, the dynamic range 
is 597 which fits (26). 

B. SPECIAL CASE: MODULI SPACING 3 OR MORE APART 

The following results have been observed, when the channel numbers are 3, 4 and 
5 distance apart as shown in Table 3.3, 3.4 and 3.5. 



nij 


m2 


Dynamic Range 


l sl Difference 


5 


8 


34 


18 


8 


11 


52 


18 


11 


14 


70 


18 


14 


17 


88 


18 


M RSN s = 3m, + 3m 2 - 5 



Table 3.3: The dynamic range of mi and m 2 = /«/ + 3 where nu > 5 (A = 3) 
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mi 


m 2 


Dynamic Range 


1 st Difference 


7 


11 


48 


24 


11 


15 


72 


24 


15 


19 


96 


24 


19 


23 


120 


24 


Mrshs = 1 + 3ni2 - 6 



Table 3.4: The dynamic range of mi and m 2 = /«/ + 4 where m / >7 (A = 4) 



mj 


m2 


Dynamic Range 


1 st Difference 


7 


12 


50 


3Q 


12 


17 


80 


30 


17 


22 


110 


30 


22 


27 


140 


30 




M rsns = 3m 1 


+ 3m 2 - 7 





Table 3.5: The dynamic range of m/ and m 2 = mj + 5 where mj > 7 (A - 5) 



These relationships show a uniform increase of 3 with each unit increase of 
value in mi and m 2 . The formulas 



M RSNS 


= 3m 1 + 3m 2 - 5 


(27) 


M RSNS 


= 3mi + 3m 2 - 6 


(28) 


MrsnS 


= 3m 1 + 3m 2 - 7 


(29) 



will work for the first, second and third sets of cases which are shown in Table 3.3, 3.4 
and 3.5. respectively.. 

As in the earlier example these cases can be expressed in terms of the lower value, 
mi and a difference, A - m 2 - ntj. In the first set A = 3, in the second set A = 4, and in the 
third set A = 5. Thus, (27), (28) and (29) can be written as 
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M rsns = 3m i + 3(mj + A) - 5 = 6m; + 4 (30) 

A 

Mrsns = 3mi + 3(m; + A) - 6 = 6m; + 6 (31) 

M r5w5 = 3m; + 3(m; + 4) - 7 = 6m; + <3. (32) 

These can be combined as 

Mrsns = 6m i + 2A- 2 (33) 

and now A can be replaced with m 2 - m;, yielding 

Mrsns ~ 6m j + 2(m 2 -m;) - 2 = 4mj + 2m? - 2. (34) 



and this formula works for all three values of A, A = 3, A = 4, and A = 5. This formula is 
checked against the computer generated results. For example, when nii=101 and 
m 2 =705, the dynamic range is 612 which fits (34). It is conjectured that (34) is valid for 
any A >3. Table 3.6 shows a summary of the specific 2-channel dynamic range formulas. 



Spacing between moduli {A ) 


Dynamic Range Formula 


A= 1 or A = 2 


Mrsns ~ 4m; + 2m 2 - 5 


A >3 


Mrsns = 4m ; + 2m 2 - 2 



Table 3.6: The Summary of the Dynamic Range Formulas for specific 2-channel 

RSNS 



C. GRAPHICAL ANALYSIS 

The algorithm rsns.m also finds and displays the point dynamic range of the 2- 
channel RSNS as explained in Chapter II. The point dynamic range versus the point 
index for the channel and shift values /«; = 4, sj = 0 and m? = 5, s? = / are shown in 
Figures 3.1 and 3.2, respectively. The waveform with the points A, B, C. D and E, that 

p 

are shown in Figures 3.1 and 3.2, construct a special waveform that repeats itself ■ RS ^ - S . 

This special waveform, the point dynamic range (PDR) waveform, gives the dynamic 
range starting at each point with increasing point index. The fundamental period of this 
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particular 2- channel RSNS example is PFrsns = 80 by (16). So, this PDR waveform 
repeats itself every 40 points. 

The maximum value of the point dynamic range is the dynamic range of the 2- 
channel RSNS. In this example, M RSNS = 21. The dynamic range point indexes are 5, 8, 
45, 48, ... for s/ = 0, s 2 = 1 and 14, 17, 54, 57, ... for si = 1 , S 2 = 0 which can be seen 
from Figure 3.1 and 3.2. These values are the beginning point index of the dynamic 
range. The dynamic range point index values are also half of the fundamental period 
apart from each other. 

The PDR waveform of the point dynamic range of the 2-channel RSNS will be the 
same for different shift values as shown in Figure 3.1 and 3.2. The only difference is the 
point index (beginning value) of the dynamic range. 



27 




Figure 3.1: The Point Dynamic Range Waveform for m t = 4, s/ = 0,and m 2 = 5, 

s 2 = 1. 



28 



The Point Dynamic Range 




Figure 3.2: The Point Dynamic Range Waveform for m } = 4, s/= 1 , and m 2 = 5, 

s 2 = 0. 
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IV. 3-CHANNEL ROBUST SYMMETRICAL NUMBER SYSTEM 



A. SPECIAL CASE: MODULI SPACING 1 APART 

In this chapter, the 3-channel RSNS is addressed. An equal distance separates the 
moduli. When the PRP channel numbers are 1 apart and the set of the 3-channel RSNS 
are increasing uniformly, the following results have been observed as shown in Table 4. 1 . 



m, 


m2 


m3 


Dynamic Range 


1 st Difference 


2 nd Difference 


3 


4 


5 


43 


39 

51 




5 


6 


1 


82 


12 • 


7 


8 


9 


133 


12 


63 


9 


10 


11 


196 










>\ 

M 

11 RSNS 


3 , 15 

= 2^ + + 7 





Table 4.1: The Dynamic Range of my, m2 = mi + 1, and m3 = mi + 2 where my >3 



The first difference (e.g., 82 — 43 = 39) and second difference (e.g., 51 - 39 = 12) 
of the dynamic range of the 3-channel RSNS for mi, m2 = mi + 1, and m3 = mi + 2 where 
mi > 3 are shown in Table 4.1. The relationship between the dynamic range differences 
can be expressed as a 2 nd degree polynomial 

p(n) = an 2 + bn + c (35) 

where n is the lowest channel value and a, b, and c are a constant real numbers. Each set 
of 3-moduli RSNS values shows a uniform increase of 2. Therefore, if the first set of 
channel values corresponds to (35), then the second set of channel values corresponds to 

p(n+2) = a(n + 2) 2 + b(n + 2) + c (36) 

= a( n 2 + 4n + 4) + bn + 2b + c 
= an 2 + 4an + 4a + bn + 2b + c 
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= p(n) + 4an + 4a + 2b. 

The first difference of the sets ((36) - (35)) is 

q(n) = p(n + 2) - p(n) = 4an + 4a + 2b. (37) 

If the first difference of the first two sets is (37), the first set difference of the second two 
sets can be written as 



q(n + 2) = 4a(n + 2) + 4a + 2b 
= 4an + 8a + 4a + 2b 
= q(n) + 8a. 

The second difference of the sets are unique and 

q(n+2) - q(n) = 8a. 

The constant “ a ” can be found by using (39) 

q(n+2) - q(n) = 8a =12 
12 3 



The constant “b" can be found by using (37) for n = 3 ( mi = 



p(n + 2) — p(n) = 39 = 4 
39 = 



U 

36+12 



i A 



3 + 4 



V2y 

+ 2b 



3) 

+ 2 b 



b 



15 

2 



The constant “c” can be found by using (35) for n = 3 (mi = 3) 

3 2 15 3 2 15 

p( 3) = 43 = — n 2 + — n + c = — 3" + — 3 + c 
2 2 2 2 

c = 1 



Finally, the closed-form expression for the dynamic range can be conjectured as 



(38) 



(39) 

(40) 



(41) 



(42) 



M 



RSNS ~ 2 fn> ^ 



15 

7 m,+7 



(43) 



where m / > 3. This formula is checked against the computer results. For example, when 
in/ = 23. m 2 = 24. and m 3 = 25, the dynamic range is 973 which fits (43). 
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B. SPECIAL CASE: MODULI SPACING 2 APART 

When the PRP channel numbers are 2 apart, the following results have been 
observed as shown in Table 4.2. 



mi 


m2 


m3 


Dynamic Range 


1 st Difference 


2 nd Difference 


5 


1 


9 


104 


126 




9 


11 


13 


230 


48 


174 


13 


15 


17 


404 


48 


222 


17 


19 


21 


626 










M RSNS Z 


3 , 21 

= 2 m > + —mi +14 





Table 4.2: The Dynamic Range of mi, m2 — mi + 2, and m3 = mi + 4 where mi > 5 



The first difference (e.g., 230 - 104 = 126) and second difference (e.g., 174 - 126 
= 48) of the dynamic range of the 3-channel RSNS for m u m2 = m; + 2, and m3 = m/ + 4 
where mi > 5 are shown in Table 4.2. The relationship between the dynamic range 
differences can also be expressed as a 2 nd degree polynomial 

p(n) = an 2 + bn + c (44) 

where n is the lower channel value and a, b, and c are constant real numbers. Each set of 
the 3-moduli RSNS values shows a uniform increase of 4. Therefore, if the first set of 
channel values corresponds to (44), then the second set of channel values corresponds to 

p(n+4) = a(n + 4) 2 + b(n + 4) + c (45) 

= a( n 2 + 8n + 16) + bn + 4b + c 
= an 2 + 8an + 16a + bn + 4b + c 
= p(n) + 8an + 16a + 4b. 

The first difference of the sets ((45) - (44)) is 

q(n) = p(n + 2) - p(n) = 8an + 16a + 4b. (46) 



33 



If the first difference of the first two sets is (46), the first set difference of the second two 
sets can be written as 



q(n + 4) = 8a(n + 4) + 16a + 4b 
= 8an + 32a +16a + 4b 
= q(n) + 32a. 

The second difference of the sets are unique and 

q(n+2) - q(n) - 32a. 

The constant “a” can be found by using (48) 

q(n+4) - q(n) = 32a =48 

48 3 

a " 32 “ 2 

The constant can be found by using (46) for n = 5 (mj = 5) 






pin + 4) - p{n) = 126 = 8 



5+16 






v2; 



+ 4 b 



120 + 48 

126 = + 4 b 

2 

b = — 

2 

The constant “c” can be found by using (44) for n = 5 (mj = 5) 

3 21 3 21 

pi 5) = 1 04 = ~n~ + ~~ n + c = ~~5~ + ~~5 + c 
2 2 2 2 

c= 14 



(47) 



(48) 

(49) 



(50) 



(51) 



Finally, the closed-form expression for the dynamic range can be conjectured as 



M 



3 2 21 

rsns ~ 2 !Ul 2 14 



(52) 



where mi >3. This formula is checked against the computer results. For example, when 
mi = 25, m 2 = 21 , and = 29, the dynamic range is 1214 which fits (52). 

C. GRAPHICAL ANALYSIS 

The algorithm (rsns.m) also finds and displays the point dynamic range of the 3- 
channel RSNS as explained in Chapter II. The point dynamic range versus the point 
index for the channel values m; = 3, m 2 = 4, and = 5 are shown in Figures 4.1 through 
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4.6. The waveform with the points A, B, C, D, E, F, G, H, I, J, K, L, and M, that are 
shown in Figures 4.1 through 4.6, construct a special waveform that repeats itself every 

p 

™ — . This special waveform, the point dynamic range (PDR) waveform, gives the 

dynamic range starting at each point. The fundamental period of this particular 3-channel 
RSNS example is 360 by (16). So, this spike waveform repeats itself every 180 points. 

The maximum value of the point dynamic range is the dynamic range of the 3- 
channel RSNS. In this example, the dynamic range is 43. The first point index that 
yields the dynamic range is 61, 142, 8, 134, 180, and 81 for Figures 4.1 through 4.6, 
respectively. These values are the beginning point index of the dynamic range. The 
dynamic range point index values are also half of the fundamental period apart from each 
other. 

The PDR waveform of the point dynamic range of the 3-channel RSNS will be the 
same for different shift values as shown in Figures 4.1 through 4.6. The beginning value 
of the dynamic range is different from each other. When one of the shift values is fixed, 
the other two shift value combinations of the PDR waveforms are opposite of each other. 
The PDR waveform of the Figures 4.1, 4.3, and 4.5 are the opposite of Figures 4.2, 4.4, 
and 4.6, respectively. 
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Figure 4.1: The Point Dynamic Range Waveform for mj = 3, si = 0, m 2 = 4, S 2 = 1, 

m 3 = 5, sj = 2. 
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The Point Dynamic Range 




Figure 4.2: The Point Dynamic Range Waveform for mj = 3, Sj = 0, m 2 = 4, S 2 = 2 , 

m 3 - 5, Sj = 1. 
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Figure 4.3: The Point Dynamic Range Waveform for nij = 3, sj = 1, m 2 = 4, s 2 = 0, 

m 3 = 5, S 3 = 2. 
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Figure 4.4: The Point Dynamic Range Waveform for nil = 3, sj = 1, m 2 = 4, s 2 = 2 , 

m 3 = 5, sj - 0. 
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The Point Dynamic Range 




m 3 = 5, s 3 = 1. 
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Figure 4.6: The Point Dynamic Range Waveform for - 3, si = 2, m 2 = 4, s 2 = 1, 

m 3 = 5, s 2 = 0. 
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V. APPLICATION OF THE RSNS TO FOLDING ANALOG-TO- 

DIGITAL CONVERTERS 



A. BACKGROUND FOR THE FOLDING ADC 

Analog to digital converters (ADCs) are critical building blocks in a wide range of 
hardware, from radar and electronic warfare systems to multimedia based personal 
computer and work stations. The best-known architecture for high-speed ADCs is the 
flash converter structure [1.5]. In this structure an array of comparators amplitude analyze 
the input voltage with a set of increasing reference voltages. The comparator outputs 
represent the input signal in a thermometer code, which can be easily converted into a 
gray or binary code. This architecture requires ( 2 N - 7) comparators to achieve TV-bit 
resolution, which makes it difficult to achieve a high resolution while maintaining, at the 
same time, a large analog bandwidth, a low power dissipation, and a small die area. 

The need constantly exists for converters with higher resolution, faster conversion 
speeds, and lower power dissipation. To reduce the number of power consuming 
components, high performance ADCs employ a parallel configuration of analog folding 
circuits to symmetrically fold the input signal prior to quantization by high speed 
comparators (analog preprocessing) [6-9]. Recently a SNS preprocessing technique has 
been described that can easily be incorporated into the established techniques to provide 
an enhanced resolution capability with a fewer comparators loaded in parallel. The 
approach is based on preprocessing the analog signal with symmetrical number systems 
(SNS and OSNS). The preprocessing is used to decompose the analog amplitude analyzer 
operation into TV sub-operations (channels) which are of smaller computational 
complexity. Each sub-operation symmetrically folds the analog signal with folding 
period equal to the channels or twice the channels. Thus, each sub-operation only 
requires a precision in accordance with that modulus. A much higher resolution is 
achieved after the TV different channels are used and the results of these low precision sub- 
operations arc combined. Increasing the resolution of these architectures can be 
accomplished by incorporating more efficient encoding. One of the important factors 
which helps to measure the efficiency of the SNS encoding is the dynamic range. The 
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maximum number of consecutive vectors without a redundancy. The other factor is the 
total number of comparators, which are required for the each encoding. The 7-bit OSNS 
ADC block diagram is shown in Figure 5.1. The problem with the SNS and OSNS is the 
encoding errors. The transfer function of 7-bit OSNS and the encoding errors for 0% 
LSB decimation width is shown in Figure 5.2. 

Folded Output 




Figure 5.1: The 7-bit OSNS Block Diagram 




Figure 5.2: The Transfer Function with Isolation Width = 0% LSB for 7-bit OSNS 
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B. IMPLEMENTATION ISSUES 



To demonstrate the efficiency of the RSNS preprocessing a unipolar 5-bit ADC 
with mi = 3, m 2 = 4, and m 3 = 5 is considered. The dynamic range of this 3-channel 
RSNS is 43. A schematic diagram of this folding ADC is shown in Figure 5.3. 

The input signal is folded in parallel with each folding period, Prsns equal to 2Nm, 
where N is the number of channels and m is the modulus (channel value). The output 
from each folding circuit is quantitized with m comparators. For this simulation, there 
will be 3 + 4 + 5= 12 comparators total and a maximum 5 of them will load in parallel. 
The comparator threshold levels are adjusted to midlevel quantize the folded output 
waveform. One interconnecting single folding and comparator stage is shown in Figure 
5.4 and the corresponding values of the resistors, current and voltage sources are shown 
in Table 5.1. The encoder to a binary representation for a high-resolution digital output 
converts the comparator outputs (the RSNS representation of the input signal). 

Each RSNS folding circuit must symmetrically fold the input signal at least 

~ 2N~ num b er t * mes - High performance folding circuits typically consists of several 

identical but independent stages inter-connected in parallel with one stage for each 
required fold. An additional stage is sometimes used to register the input voltages in DC. 

V re f is a supplied reference voltage to each stage equal to the input voltage point at 
which a fold (peak) is desired. Also shown are the analog input node and the folded 
output node. For this 5-bit example, a LSB code width of 12mv with a full-scale voltage 
(2 6 -l) x 12mv = 0.756 v is chosen. The folding period for the modulus channel is 



to N. The first required reference voltages (V refur above zero and V re f(or below zero) are 



T m =2Nm i x(LSB Code Width ) v 



(53) 



where N is the number of modulus and ntj is the value of the modulus and it goes from 1 




( 54 ) 
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Figure 5.3: The 5-bit RSNS ADC Block Diagram 
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Stage 1 




Figure 5.4: One Interconnecting Single Folding and Comparator Stages of the 5-bit 

RSNS 
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Resistor Design 


R3 


400 0 


Values 


R4 


300 0 




R5 


300 0 




R6 


250 0 




R7 


250 O ( mod-3) - 300 O (mod-4) - 350 O (mod-5) 


Voltage Sources 


VI 


15 v 


Design Values 


V2 


15 v 




V3 


15 v 


Current Sources 


11 


0.001 A 


Design Values 


12 


0.001 A 



Table 5.1: The Design Values of the 5-bit RSNS (ni/ = 3, m 2 = 4, m 3 =5) 



T 

V„ /(0) = — ^-+ ( s t x(LSB Code Width)) v (55) 
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where T m is the folding period and s, is the shift values of the modulus nij. The other 

required reference voltages are 

T 

V ref(N) = -^+ (j,x (LSB Code Width)) + NT mj v (56) 

where N is the N 11 ' reference voltages and lies between {0, 1,2, . . . } for the modulus m ( . 

The folding period by (53) for the mi = 3 is 7j = 2 x 3 x 3 x 12mv = 216 mv. That 
is, the folding period is scaled to the modulus. For mi = 5, the shift value is taken si = 0. 

T i 

The first required reference voltage, above zero is then V refm = — = 108 mv by (54). The 

T 3 

reference voltage, below zero is V ref(0) = — ^- = -108mv by (55). The other required 

reference voltages for mj = 3 are calculated by (56) and shown in Table 5.2. 

The folding period by (53) for the m 2 = 4 is T 4 = 2 x 3 x 4 x 12mv = 288 mv. That 
is, the folding period is scaled to the modulus. For m 2 = 4, the shift value is taken S 2 = 1. 
The first required reference voltage, above zero is then 

T 

V rem =-^- + (s l x(LSB Code Width)) = \56mv by (54). The reference voltage, below 

T 

zero is V ref(0) = — -^- + {s 2 X{LSB Code Width)) = -132/«v by (55). The other required 

reference voltages for m 2 = 4 are calculated by (56) and shown in Table 5.2. 

The folding period for the m 3 = 5 is T 5 = 2 x 3 x 5 x 12mv = 360 mv by (53). That 
is, the folding period is scaled to the modulus. For mj = 5, the shift value is taken S 3 = 2. 
The first required reference voltage, above zero is 

T 

then V ref(l) =-^- + (s 3 x(LSB Code Width)) = 204mv by (54). The reference voltage, 

T 

below zero is V ref(0) = — -^- + (s } x(LSB Code Width)) =-156mv . The other required 
reference voltages for m 3 = 5 are calculated by (56) and shown in Table 5.2. 
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Fold 


Modulus 3 


Modulus 4 


Modulus 5 


Number 


(' T 3 =0.216 ) 


(7> 0.288) 


(7>= 0.360) 


0 


-0.108 


-0.132 


-.0156 


1 


0.108 


0.156 


0.204 


2 


0.324 


0.444 


0.564 


3 


0.540 


0.732 


0.924 


4 


0.756 


1.020 


1.284 


5 


0.972 


1.308 


1.644 


6 


1.188 


1.596 


- 


7 


1.404 


- 


- 


8 


1.620 







Table 5.2: The Reference Voltages of the 5-bit RSNS Design 



A SPICE simulation is used to generate the folding waveforms and comparator 
outputs states. A logic map (or an algorithm (spice.m) to give a decimal output as shown 
in the Appendix processes the sampled folding waveforms and comparator states. The 
steady-state folding waveforms for each channel are shown in Figures 5.5, 5.6, and 5.7. 
These folding waveforms are generated using a step size of Av = 0.001 v in the SPICE 
simulation. 
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The comparator threshold voltages V, for each modulus derived from the folding 
waveforms to mid-level quantize the input signal into RSNS format. The threshold levels 
occur at the code transition points and are tailored to the particular folding waveform 
being instrumented. The comparator threshold voltages for the 5-bit RSNS are shown in 
Table 5.3. Note the threshold values are not uniformly spaced. The comparator outputs 
are shown in Figures 5.5, 5.6 and 5.7 with their corresponding folded waveform for 5-bit 
RSNS. 



Comparator 

Number 


Modulus-3 


Modulus-4 


Modulus-5 


1 


5.4097 


5.2994 


5.2234 


2 


5.6748 


5.5169 


5.4045 


3 


5.9452 


5.7408 


5.5946 


4 


- 


5.9668 


5.7878 


5 


- 


- 


5.9818 



Table 5.3: The Comparator Threshold Voltages for the 5-bit RSNS Design 
C. TRANSFER FUNCTION 

To obtain a more convenient digital output, the RSNS encoded signals from each 
channel are recombined. The comparator outputs from each channel are decoded using a 
logic block or programmable logic array (PLA). The RSNS-to-decimal mapping function 
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Figure 5.5: Folding Output Waveform/Comparator Outputs for the RSNS m; = 3 



51 



MOD4 




0.3 0.4 0.5 0.6 0.7 0.8 



20 1 1 1 1 1 1 1 r 




_20 1 1 1 1 i i i i 1 1 1 

-0.2 -0.1 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 



20 1 1 1 1 1 1 1 1 1 r 



CM n 

O o- 



- 20 1 1 1 1 1 1 1 1 1 1 1 

-0.2 -0.1 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 




O 



-20 

0.2 0.3 0.4 

ADC INPUT VOLTAGE 

Figure 5.6: Folding Output Waveform/Comparator Outputs for the RSNS m 2 = 4. 
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•e 5.7: Folding Output Waveform/Comparator Outputs for the RSNS mj = 5 . 
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for the 5-bit RSNS design is shown in Table 5.4 where each entry indicates the number of 
comparators in the ON state. 



The Point Index 


RSNS Moduli 


3 


4 


5 


0 


0 


0 


0 


1 


1 


0 


0 


2 


1 


1 


0 


3 


1 


1 


1 


4 


2 


1 


1 


5 


2 


2 


1 


6 


2 


2 


2 


7 


3 


2 


2 


8 


3 


3 


2 


9 


3 


3 


3 


10 


2 


3 


3 


11 


2 


4 


3 


12 


2 


4 


4 


13 


1 


4 


4 


14 


1 


3 


4 


15 


1 


3 


5 


16 


0 


3 


5 


17 


0 


2 


5 


18 


0 


2 


4 


19 


1 


2 


4 


20 


1 


i 


4 



Table 5.4: The 5-bit RSNS-to-Decimal Mapping 
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The RSNS ADC architecture is simulated using a sampling period of Av = 0.001 v 
with the corresponding transfer function shown in Figure 5.8. The LSB code width is 
assumed 12 mv for this simulation. Therefore, the dynamic range lies between the input 
voltages of 0.012 x 61 = 0.732 v and 0.012 x 103 = 1.236 v. The entire graph. Figure 5.2 
lies in the dynamic range. This corresponds to the portion of Figure 5.8 where the 
dynamic range lies from 0.732 v to 1.236 v. When the sampling period is decreased to Av 
= 0.0004 v, the transfer function does not change. Note that the encoding errors typically 
present in the SNS and OSNS ADCs do not exist. This clearly demonstrates the 
robustness of the RSNS encoding (Gray code properties). For this system, channel and 
shift values, the dynamic range lies between the point indexes 61 and 103 as shown in 
Table 2.4. 

The comparator states of this architecture are simulated to the output decimal 
value by using the algorithm, trans.m. The beginning point index (61) of the dynamic 
range is assigned the decimal value 0. The dynamic range point indexes are assigned the 
decimal values from 0 to 42. The point indexes within the dynamic range are assigned 
the decimal values starting from 0 to 42. Then, starting from the point index 60 to 0, the 
decimal values are assigned backwards. If the point index of the discrete states is a 
redundancy within the dynamic range, it is assigned the corresponding decimal value. If 
not, it receives a new decimal value starting from -1. Finally, the same procedure is 
applied to the point index starting from 104, but the values of the decimal codes start 
from 43. The dynamic range (useful code combinations) and the redundancies 
(ambiguities) can be seen from Figure 5.8. 
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VI. CONCLUDING REMARKS 



The comparison of the RSNS with the other SNS formulations is shown in Table 
6. 1 . The folding period for the RSNS is a factor of 2N greater then the SNS and N greater 







Number of 




Scheme 


Folding Periods 


Comparators 
Per Channel 


Dynamic Range 


SNS 


m 




m, 

. 2 . 




N 
i = 1 


SNS 








(all moduli 
odd) 


m 




m, 

. 2 . 






SNS 








(one moduli 
even) 


mi 




>», 

_ 2 _ 




+ r i m i\ 

[ * 1=2 / = >+ 1 J 


OSNS 


2nij 


mi - . 


I 


N 

tl m 

/= 1 


RSNS 






4mi + 5m? — 5 


(2 moduli) 


2Nm, 


mi 


(moduli 1 or 2 apart) 


RSNS 






4m/ + 5ni2 - 2 


(2 moduli) 


2Nnii 


m 


(moduli 3 or more apart) 


RSNS 
(3 moduli) 


2Nnii 


m { 


3 2 15 

j'". +Y m i + 7 








(moduli 1 apart) 


RSNS 
(3 moduli) 


2Nm, 


m t 


3 , 21 

—nt+ — nt + 14 

2 -i 2 i 








(moduli 2 apart) 



Table 6.1: The Summary of the SNS Systems 
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then the OSNS. The number of comparators per channel for the RSNS is the maximum of 
all the SNS systems. The dynamic range for the RSNS is less than the OSNS, but greater 
than the SNS. 

The essential contribution of this thesis is defining a new symmetrical number 
system, robust symmetrical number system (RSNS) and investigating its properties. High 
Speed Analog-to-Digital Converters (ADCs) for 2-channel and 3-channel are summarized 
in Tables 6.2 and 6.3, respectively. 



Number 

System 


Channel 

(Modulus) 


Dynamic 

Range 


Folding 

Period 


Number of 
Comparators 


SNS 


4 


1 

S' 


4 


2 


5 


5 


. 2 


OSNS 


4 


20 


8 


3 


5 


10 


4 


RSNS 


4 


21 


16 


4 


5 


20 


5 



Table 6.2: 2-Channel (modulus) Summary of the Number Systems for m.i-4 and 

m 2 = 5 

The dynamic range of the 2-channel RSNS is greater than the SNS and OSNS, but 
it has an increased number of comparators. For the 3-channel RSNS, although the 
dynamic range is greater than the SNS, it is less than the OSNS. The RSNS has Gray code 
properties in its behavior. Gray code properties mean only one comparator threshold is 
crossed at a time across parallel channels (moduli). This property of the RSNS eliminates 
the possibility of large errors at the code transition points and the requirements for 
interpolation signal processing. 

Using SPICE, explained in Chapter V, simulates the 3-channel RSNS for specific 
channel and shift values. From the simulation results, it can be seen that the RSNS does 
not have any encoding errors. This property of the RSNS makes itself an attractive and 
useful number system in High Speed ADCs. 
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Number 

System 


Channel 

(Modulus) 


Dynamic 

Range 


Folding 

Period 


Number of 
Comparators 


SNS 


3 


11 


3 


1 


4 


4 


2 


5 


5 


2 


OSNS 


3 


60 


6 


2 


4 


8 


3 


5 


10 


4 


RSNS 


3 


43 


18 


3 


4 


24 


4 


5 


30 


5 



Table 6.3: 3-Channel (modulus) Summary of the Number Systems for mi = 3, m 2 = 4 

and mi = 5 

Future efforts will attempt to resolve the following important questions. 

• How does the shift effect the beginning point index of the dynamic range for 
any RSNS configuration? 

• Can a general closed-form expression be developed for the 3-channel and N- 
channel RSNS? 
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APPENDIX 



A . DESCRIPTION OF THE COMPUTER ALGORITHM DEVELOPED 

The programming in this study is accomplished using the MATLAB software 
package. The main program, which is used to find the properties of the RSNS, is the 
rsns.m. The program spice.m and trans.m are used to plot the SPICE simulation. 

The program is rsns.m finds and displays the dynamic range, fundamental period 
and beginning-ending point index of the discrete states without redundancy of any 
number of RSNS channels. The channel and shift values must be entered to run the 
program. This program also finds and displays the point dynamic range of the RSNS. 

The program spice.m reads the folding waveform and comparator outputs data and 
displays the folding waveform and comparator outputs for each channel (modulus) of the 
RSNS by using SPICE simulation. 

The program trans.m also reads the folding waveform and comparator output data 
which are calculated by using SPICE simulation. The algorithm, trans.m finds the 
decimal output codes for each comparator outputs of the RSNS and displays the transfer 
function of it. The dynamic range and its location can be seen in this algorithm display 
(figure). 
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B. PROGRAM LISTING 



1. RSNS Properties (rsns.m) 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 



% rsns.m % 

% Modified July 05 th 1 997 % 

% This program finds the dynamic range, point dynamic range, fundamental period % 
% and beginning-ending values of ambiguity vectors for N-channel RSNS % 



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

clear all 

% Enter the channel number and the values of them 

disp('This program finds the maximal strings of non-redundant vectors'); 

disp('for the N channel ROBUST SNS'); 

% Define a variable for the channel number 

chanum=input('Enter the Number of Channels for ROBUST SNS » '); 

% Find the fundamental period with respect to the entered channel and shift values 

period= 1 ; 

for i=l:chanum 

m(i)=input(['Enter ' int2str(i) '.Channel Value » ']); 
period=LCM(period,(2*chanum*m(i))); 
end 

nsearch=period+30; 

prompt='y'; 

% Enter the corresponding shift values 
while (prompt=='y')l(prompt=='Y') 
for i=l:chanum 

s(i)=input( ['Enter ' int2str(i) '. Channel Shift Value » ']); 
end 
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% Initiate the variables to zero 



i=0; 

ii=0; 

j=°; 

jj=o; 

k=0; 

% Define a variable for the fundamental period 
funper=0; 

% Define a variable for the dynamic range 
dynrange=0; 

% Define the waveform for the entered channel and shift values and put the data in 
% matrix “g” 
for r=l:chanum 

% Define the matrix “mm” to put the entered channel and shift values 
mm(r,[l 2])=[m(r) s(r)]; 

% Use the definition of RSNS 
for i=l+s(r):chanum*m(r)+s(r) 
g(r,i)=floor((i-s(r))/chanum); 
end 

for i=chanum*m(r)+l+s(r): 2*chanum*m(r)+s(r) 

g(r,i)=floor((2*chanum*m(r)+chanum-i+s(r)- 1 )/chanum); 
end 

% Extent enough the data of the folding waveforms of the RSNS to find the dynamic 
% range and fundamental period 

g(r,2*chanum*m(r)+s(r)+l:4*chanum*m(r)+s(r))=g(r,l+s(r):2*chanum*m(r)+s(r)); 

g(r,4*chanum*m(r)+s(r)+l:8*chanum*m(r)+s(r))=g(r,l+s(r):4*chanum*m(r)+s(r)); 

g(r,8*chanum*m(r)+s(r)+l:16*chanum*m(r)+s(r))=g(r,l+s(r):8*chanum*m(r)+s(r)); 
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g(r, 1 6*chanum*m(r)+s(r)+ 1 :32*chanum*m(r)+s(r))=g(r, 1 +s(r): 1 6*chanum*m(r)+s(r)); 
g(r,32*chanum*m(r)+s(r)+l :64*chanum*m(r)+s(r))=g(r,l+s(r):32*chanum*m(r)+s(r)); 
g(r,64*chanum*m(r)+s(r)+l: 128*chanum*m(r)+s(r))=g(r,l+s(r):64*chanum*m(r)+s(r)); 
end 

% Define a new matrix “ga” which is the transpose of matrix “g” 
ga=g’; 



% Define a new matrix “gb” which is extended one more column to show the number of 
% rows 

gb(:,[2:(chanum+l)])=ga(:,[l:chanum]); 

[sgbr,sgbc]=size(gb); 

gb(:,l)=(l:l:sgbr)’; 

% Define a new matrix “gc” that gives the row number through the number of search 
gc=gb; 

gc(nsearch+l:sgbr, :)=[]; 

[sgcr,sgcc]=size(gc); 

% Find the fundamental period to check this result with the formula result of it 
xper=gc(l,[2:(chanum+l)]); 
for jj=chanum:sgcr; 

if gc(jj,[2:(chanum+l)])==xper 
funper=gc((jj-(chanum- 1 )), 1 ); 
break 
end 

end 

% Find the first redundancies in matrix “gc” 
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k=l; 

for ii=chanum- 1 insearch; 
xrec=gc(ii , [2 : (chanum+ 1 )] ); 
for ij=l :floor((nsearch-ii)/(2*chanum)) 
jj=ii+2*chanum*ij; 

if gc(jj,[2:(chanum+l)])==xrec 
redun=gc(jj,l); 

% Define a new matrix “h” to put the first redundancies 
h(k,l)=ii; 
h(k,2)=redun; 
k=k+l; 
break 
end 

end 

end 

% Define a new matrix “ hsort” to sort the redundancy matrix “h” 
hsort=h; 

[yoy,ioi]=sort(hsort); 

% Define a new matrix “ hsorted” to put the actual row number of “hsort” 

hsorted=[yoy(ioi(:,2), 1 ) yoy(:,2)]; 

hsorted; 

% Define a new matrix “hreduced” to eliminate the rows of the matrix “hsorted” and that 

% will not allow the first column to be monotone increasing 

[ssr,ssc]=size(hsorted); 

hreduced=hsorted; 

a=h(ssr,l); 
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[rx cx]=size(hreduced); 
for k=l:ssr 
for i=l:ssr 
if i<rx 

if hreduced(i, 1 )==a 
hreduced(i+l :ssr,:)=[]; 
break 

elseif hreduced(i+ 1 , 1 )<hreduced(i, 1 ) 

hreduced(i+l, :)=[]; 

break 

end 

end 

[rx cx]=size(hreduced); 
end 
end 

hreduced; 

% Define a new matrix “H” that shows the point index of the distinct vectors without 
% redundancy and their lengths 
[hsr,hsc]=size(hreduced); 

H(l,l)=(chanum-1); 

H(2:hsr+ 1 , 1 )=hreduced( 1 :hsr, 1 )+ 1 ; 

H( 1 :hsr,2)=hreduced( 1 :hsr,2)- 1 ; 

H(hsr+ 1 ,2)=nsearch; 

H( 1 :hsr+ 1 ,3)=H( 1 :hsr+ 1 ,2)-H( 1 :hsr+ 1 , 1 )+ 1 ; 

[Hsr,Hsc]=size(H); 

% Find the point dynamic range of the RSNS and put these values to the matrix “HH” 
t=0; 
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for n=l:Hsr 

for p=H(n,l):H(n+l,l) 

HH(p, 1 )=p; 

HH(p,2)=H(n,2); 

HH(p,3)=H(n,3)-t; 

t = t+ 1 

if t==(H(n+ 1 , 1 )-H(n, 1 )+ 1 ) 
t=0; 
end 
end 

n=n+l; 
if n==Hsr 
break 
end 
end 

[HHsr,HHsc]=size(HH); 

HH(HHsr,:)=[]; 

% Plot the point dynamic range for the entered channel and shift values of RSNS 
figure(l) 

plot(HH(:,l),HH(:,3)),grid 

xlabel('THE SHIFTED CHANNEL VALUES’),ylabel(’THE DYNAMIC RANGE OF 
THE CHANNEL VALUES’) 

title(['THE DYNAMIC RANGE DISTRIBUTION OF,num2str(chanum),'-CHANNEL 

RSNS']) 

orient tall 

%print 

% Find the dynamic range of the entered channel and shift values of RSNS 
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HHH=max(H); 

dynrange=HHH(:,3); 

% Display the dynamic range, fundamental period, entered channel and shift values and 
% matrix “H” which shows the % beginning-ending point index of the distinct vectors 
% without redundancy and their lengths 
dispC ') 

disptf’THE FUNDEMANTAL PERIOD IS ',num2str(funper),' ’]) 
dispC ') 

disp(['THE DYNAMIC RANGE IS ',num2str(dynrange),' ']) 

fprintfCVn BEGIN POINT INDEX END POINT INDEX DYNAMIC RANGE\n') 

fprintf(’% 1 1 .Of % 1 6.0f % 1 2.0f \n',H’) 

fprintf('\n CHANNEL VALUES SHIFT VALUES\n') 

fprintf(’% 1 1 .Of % 1 7.0f Vn’.mm’) 

prompt=input('Would you like to try another shift (y/n) ? »',’s'); 
end 

2. Digital Signal Processing for the RSNS ADC (spice.m) 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 



% spice.m % 

% Modified July 10 th 1997 % 

% This program finds the folded waveform and comparator outputs of % 

% the 3-channel RSNSwhich is simulated in SPICE program % 

% mj = 3 si = 0 % 

% m 2 -4 S 2 = J % 

% mj = 5 sj = 2 % 



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

clear all 

% Load the folding circuit and comparator outputs of the modulus-3 RSNS 
load newmod3.dat 
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a=newmod3; 



% Plot the folding waveform and comparator outputs for the modulus-3 RSNS 
figure(l) 

subplot(4,l,l),plot(a(:,l),a(:,2),'y'),grid,title('MOD3'),ylabel('FOLDED OUTPUT’) 
subplot(4, 1 ,2),plot(a(: , 1 ),a( : ,5),'r'),grid,ylabel('C3') 
subplot(4, 1 ,3),plot(a(:, 1 ),a(:,4)/g’),grid,ylabel('C2') 

subplot(4,l,4),plot(a(:,l),a(:,3),'b'),grid,ylabel('Cl'),xlabel('ADC INPUT VOLTAGE') 
orient tall 

% Load the folding circuit and comparator outputs of the modulus-4 RSNS 

load newmod4.dat 

b=newmod4; 

% Plot the folding waveform and comparator outputs for the modulus-4 RSNS 
figure(2) 

subplot(5,l,l),plot(b(:,l),b(:,2),'y'),grid,title('MOD4'),ylabel('FOLDED OUTPUT') 
subplot( 5 , 1 ,2) ,plot(b(: , 1 ) ,b( : ,6) ,'m') ,grid,ylabel('C4') 
subplot(5,l,3),p* ot (b(:,l),b(:,5),'r'),grid,ylabel('C3') 
subplot(5,L4),plot(b(:,l),b(:,4),'g'),grid,ylabel('C2') 

subplot(5,l,5),plot(b(:,l),b(:,3),'b'),grid,ylabel('Cl'),xlabel(’ADC INPUT VOLTAGE') 
orient tall 

% Load the folding circuit and comparator outputs of the modulus-5 RSNS 

load newmod5.dat 

c=newmod5; 

% Plot the folding waveform and comparator outputs for the modulus-5 RSNS 
figure(3) 
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subplot(6, 1 , 1 ),plot(c(:,l ),c(:,2)/y'),grid,title(’MOD5'),ylabel('FOLDED OUTPUT') 
subplot(6, 1 ,2),plot(c(: , 1 ),c(: ,7),'m'),grid,ylabel('C5’) 
subplot(6,l,3),plot(c(:,l),c(.\6),'r'),grid,ylabel('C4') 
subplot(6,l,4),plot(c(:,l),c(:,5),'g'),grid,ylabel('C3') 
subplot(6, 1 ,5),plot(c(: , 1 ),c(:,4), , c’),grid,ylabel( , C2 , ) 

subplot(6, 1 ,6),plot(c(: , 1 ),c(: ,3),'b'),grid,ylabel('C 1 '),xlabel('ADC INPUT VOLTAGE') 
orient tall 

% Plot the folding waveforms for modulus-3, 4, and 5 RSNS together 
figure(4) 

subplot(3, 1 , 1 ),plot(a(:, 1 ),a(:,2),'y' ),grid,title('MOD3'),ylabel('FOLDED OUTPUT') 
subplot(3,l,2),plot(b(:,l),b(:,2),'y'),grid,title('MOD4'),ylabel(’FOLDED OUTPUT') 
subplot(3, 1 ,3),plot(c(:,l ),c(:,2),'y' ),grid,title('MOD5'),ylabel('FOLDED OUTPUT 1 ) 
orient tall 

3. Transfer Function for the RSNS ADC (trans.m) 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 



% trans.m % 

% Modified July 10 th 1997 % 

% This program finds transfer function of the 3-channel RSNS, which is simulated % 
% in SPICE program. % 

% ni] = 3 S] = 0 % 

% m 2 = 4 S 2 = 1 % 

% m 3 -5 sj - 2 % 



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

clear all 
orient tall 

% Initiate the variables to zero 

i=0; 

j=0; 
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1 = 0 ; 
m=0; 
n=0; 
ka= 1 ; 
kb=l; 
kc=l; 
kd=l; 
ke= 1 ; 



% Load the folding circuit and comparator outputs of the modulus-3 RSNS and find 
% howmany comparators are on 
load newmod3.dat 
a=newmod3; 

[sar,sac]=size(a); 
for i= 1 tsar; 

al=a(i,3); 
a2=a(i,4); 
a3=a(i,5); 
if ( a 1 > 0); 

A(ka,l)=l; 

else 

A(ka. 1)=0; 
end 

if ( a2 >0); 

A(ka,2)=l; 

else 

A(ka.2)=0; 

end 

if ( a3 > 0); 
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A(ka,3)=l; 

else 

A(ka,3)=0; 

end 

ka=ka+ 1 ; 
end 
end 

% Load the folding circuit and comparator outputs of the modulus-4 RSNS and find 
% howmany comparators are on 
load newmod4.dat 
b=newmod4; 

[sbr,sbc]=size(b); 
for j=l:sbr; 

bl=b(j,3); 
b2=b(j,4); 
b3=b(j,5); 
b4=b(j,6); 
if ( bl> 0); 

B(kb,l)=l; 

else 

B(kb,l)=0; 

end 

if ( b2 > 0); 

B(kb,2)=l; 

else 

B(kb,2)=0; 

end 

if ( b3> 0); 
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B(kb,3)=l; 

else 

B(kb,3)=0; 

end 

if ( b4> 0); 

B(kb,4)= 1 ; 
else 

B(kb,4)=0; 

end 

kb=kb+l; 

end 

end 

% Load the folding circuit and comparator outputs of the modulus-5 RSNS and find 
% howmany comparators are on 
load newmod5.dat 
c=newmod5; 

[scr,scc]=size(c); 
for l=l:scr; 

cl=c(l,3); 
c2=c(l,4); 
c3=c(l,5); 
c4=c(l,6); 
c5=c(l,7); 
if ( cl> 0); 

C(kc,l)=l ; 
else 

C(kc,l)=0; 

end 
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if ( c2 > 0); 

C(kc,2)=l; 

else 

C(kc,2)=0; 

end 

if ( c3> 0); 

C(kc,3)=l; 

else 

C(kc,3)=0; 

end 

if ( c4> 0); 

C(kc,4)=l; 

else 

C(kc,4)=0; 

end 

if ( c5> 0); 

C(kc,5)=l; 

else 

C(kc,5)=0; 

end 

kc=kc+l; 

end 

end 

% Define a new matrix “D” that has every data of three folding circuit and shows how 
% many comparators are on or off in each of them 
D(:,l)=a(:,l); 
for m=l:sar; 

D(kd,2)=A(kd, 1 )+A(kd,2)+A(kd,3); 
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D(kd,3)=B(kd, 1 )+B(kd,2)+B(kd,3)+B(kd,4); 
D(kd,4)=C(kd,l)+C(kd,2)+C(kd,3)+C(kd,4)+C(kd,5); 
kd=kd+ 1 ; 
end 

% Define a new matrix “DEC” that finds the transfer function 
for ke=l:sar; 
aa=D(ke,2); 
bb=D(ke,3); 
cc=D(ke,4); 
dd=D(ke, 1 ); 

if((aa==0)&(bb==0)&(cc== 1 )) 

DEC(ke,l)=dd; 

DEC(ke,2)=-27; 

elseif((aa==0)&(bb==0)&(cc==0)) 

DEC(ke,l)=dd; 

DEC(ke,2)=-26; 

elseif ((aa==l)&(bb==0)&(cc=0)) 

DEC(ke,l)=dd; 

DEC(ke,2)=-25; 

elseif ((aa== 1 )&(bb= 1 )&(cc=0)) 

DEC(ke,l)=dd; 

DEC(ke,2)=33; 

elseif ((aa== 1 )&(bb==l )&(cc==l )) 

DEC(ke,l)=dd; 

DEC(ke,2)=34; 

elseif ((aa==2)&(bb== 1 )&(cc== 1 )) 

DEC(ke, 1 )=dd; 

DEC(ke,2)=35; 
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elseif ((aa==2)&(bb==2)&(cc== 1 )) 
DEC(ke,l)=dd; 

DEC(ke,2)=6; 

elseif ((aa==2)&(bb==2)&(cc==2)) 
DEC(ke,l)=dd; 

DEC(ke,2)=7; 

elseif ((aa==3)&(bb==2)&(cc==2)) 
DEC(ke,l)=dd; 

DEC(ke,2)=-24; 

elseif ((aa==3)&(bb==3)&(cc=2)) 
DEC(ke,l)=dd; 

DEC(ke,2)=-23; 

elseif ((aa==3)&(bb==3)&(cc==3)) 
DEC(ke,l)=dd; 

DEC(ke,2)=22; 

elseif ((aa==2)&(bb==3)&(cc=3)) 
DEC(ke,l)=dd; 

DEC(ke,2)=23; 

elseif ((aa==2)&(bb==4)&(cc==3)) 
DEC(ke,l)=dd; 

DEC(ke,2)=24; 

elseif ((aa==2)&(bb==4)&(cc==4)) 
DEC(ke,l)=dd; 

DEC(ke,2)=-22; 

elseif ((aa== 1 )&(bb==4)&(cc==4)) 
DEC(ke,l)=dd; 

DEC(ke,2)=-21; 

elseif ((aa==l)&(bb==3)&(cc==4)) 
DEC(ke,l)=dd; 
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DEC(ke,2)=-20; 

elseif ((aa==l)&(bb==3)&(cc==5)) 
DEC(ke,l)=dd; 

DEC(ke,2)=-19; 

elseif ((aa==0)&(bb==3)&(cc==5)) 
DEC(ke,l)=dd; 

DEC(ke,2)=-18; 

elseif ((aa==0)&(bb==2)&(cc==5)) 
DEC(ke,l)=dd; 

DEC(ke,2)=-17; 

elseif ((aa==0)&(bb==2)&(cc=4)) 
DEC(ke,l)=dd; 

DEC(ke,2)=-16; 

elseif ((aa==l )&(bb==2)&(cc==4)) 
DEC(ke,l)=dd; 

DEC(ke,2)=-15; 

elseif ((aa=l )&(bb==l )&(cc==4)) 
DEC(ke, 1 )=dd; 

DEC(ke,2)=15; 

elseif ((aa== 1 )&(bb== 1 )&(cc==3)) 
DEC(ke,l)=dd; 

DEC(ke.2)=10; 

elseif ((aa==2)&(bb==l)&(cc==3)) 
DEC(ke.l)=dd; 

DEC(ke,2)=41 ; 

elseif ((aa==2)&(bb==0)&(cc==3)) 
DEC(ke,l )=dd; 

DEC(ke,2)=-14; 

elseif ((aa==2)&(bb==0)&(cc==2)) 
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DEC(ke,l)=dd; 

DEC(ke,2)=37; 

elseif ((aa==3)&(bb==0)&(cc==2)) 
DEC(ke,l)=dd; 

DEC(ke,2)=38; 

elseif ((aa==3)&(bb== 1 )&(cc==2)) 
DEC(ke,l)=dd; 

DEC(ke,2)=39; 

elseif ((aa==3 )&(bb== 1 )&(cc== 1 )) 
DEC(ke,l)=dd; 

DEC(ke,2)=-13; 

elseif ((aa==2)&(bb==l)&(cc==l)) 
DEC(ke,l)=dd; 

DEC(ke,2)=35; 

elseif ((aa==2)&(bb==2)&(cc== 1 )) 
DEC(ke,l)=dd; 

DEC(ke,2)=6; 

elseif ((aa==2)&(bb==2)&(cc==0)) 
DEC(ke,l)=dd; 

DEC(ke,2)=-12; 

elseif ((aa==l)&(bb==2)&(cc==0)) 
DEC(ke,l)=dd; 

DEC(ke,2)=32; 

elseif ((aa==l )&(bb==3)&(cc==0)) 
DEC(ke,l)=dd; 

DEC(ke,2)=-l 1; 

elseif ((aa==l)&(bb==3)&(cc==l)) 
DEC(ke,l)=dd; 

DEC(ke,2)=28; 
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elseif ((aa==0)&(bb==3)&(cc== 1 )) 
DEC(ke,l)=dd; 

DEC(ke,2)=29; 

elseif ((aa==0)&(bb==4)&(cc==l)) 
DEC(ke,l)=dd; 

DEC(ke,2)=-10; 

elseif ((aa==0) &(bb==4) &(cc==2)) 
DEC(ke,l)=dd; 

DEC(ke,2)=-9; 

elseif ((aa== 1 )&(bb==4)&(cc==2)) 
DEC(ke,l)=dd; 

DEC(ke,2)=26; 

elseif ((aa== 1 )&(bb==3)&(cc==2)) 
DEC(ke,l)=dd; 

DEC(ke,2)=27; 

elseif ((aa== 1 )&(bb==3)&(cc==3)) 
DEC(ke,l)=dd; 

DEC(ke,2)=-8; 

elseif ((aa==2)&(bb==3)&(cc==3)) 
DEC(ke,l)=dd; 

DEC(ke,2)=23; 

elseif ((aa==2)&(bb==2)&(cc==3)) 
DEC(ke, 1 )=dd; 

DEC(ke,2)=42; 

elseif ((aa==2)&(bb==2)&(cc==4)) 
DEC(ke,l)=dd; 

DEC(ke,2)=19; 

elseif ((aa==3)&(bb==2)&(cc==4)) 
DEC(ke,l)=dd; 
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DEC(ke,2)=20; 

elseif ((aa==3)&(bb==l)&(cc==4)) 
DEC(ke,l)=dd; 

DEC(ke,2)=-7; 

elseif ((aa==3)&(bb==l)&(cc==5)) 
DEC(ke,l)=dd; 

DEC(ke,2)=-6; 

elseif ((aa==2)&(bb== 1 )&(cc==5)) 
DEC(ke,l)=dd; 

DEC(ke,2)=17; 

elseif ((aa==2)&(bb==0)&(cc==5)) 
DEC(ke,l)=dd; 

DEC(ke,2)=-5; 

elseif ((aa==2)&(bb==0)&(cc==4)) 
DEC(ke,l)=dd; 

DEC(ke,2)=-4; 

elseif ((aa== 1 )&(bb==0)&(cc==4)) 
DEC(ke,l)=dd; 

DEC(ke,2)=14; 

elseif ((aa== 1 )&(bb== 1 )&(cc==4)) 
DEC(ke,l)=dd; 

DEC(ke,2)=15; 

elseif ((aa==l)&(bb==l)&(cc==3)) 
DEC(ke,l)=dd; 

DEC(ke,2)=10; 

elseif ((aa==0)&(bb==l)&(cc==3)) 
DEC(ke,l)=dd; 

DEC(ke,2)=l 1; 

elseif ((aa==0)&(bb==2)&(cc==3)) 
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DEC(ke,l)=dd; 

DEC(ke,2)=-3; 

elseif ((aa==0)&(bb==2)&(cc==2)) 
DEC(ke,l)=dd; 

DEC(ke,2)=-2; 

elseif ((aa== 1 )&(bb==2)&(cc==2)) 
DEC(ke,l)=dd; 

DEC(ke,2)=8; 

elseif ((aa== 1 )&(bb==3)&(cc==2)) 
DEC(ke,l)=dd; 

DEC(ke,2)=27; 

elseif ((aa==l)&(bb==3)&(cc==l)) 
DEC(ke,l)=dd; 

DEC(ke,2)=-l; 

elseif ((aa==2)&(bb==3)&(cc== 1)) 
DEC(ke, 1 )=dd; 

DEC(ke,2)=5; 

elseif ((aa==2)&(bb==4)&(cc==l)) 
DEC(ke,l)=dd; 

DEC(ke,2)=0; 

elseif ((aa==2)&(bb==4)&(cc==0)) 
DEC(ke,l)=dd; 

DEC(ke,2)=l; 

elseif ((aa==3)&(bb==4)&(cc==0)) 
DEC(ke,l)=dd; 

DEC(ke,2)=2; 

elseif ((aa==3)&(bb==3)&(cc==0)) 
DEC(ke,l)=dd; 

DEC(ke,2)=3; 
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elseif ((aa==3)&(bb==3)&(cc== 1 )) 
DEC(ke,l)=dd; 

DEC(ke,2)=4; 

elseif ((aa==2)&(bb==3)&(cc== 1 )) 
DEC(ke, 1 )=dd; 

DEC(ke,2)=5; 

elseif ((aa==2)&(bb==2)&(cc== 1 )) 
DEC(ke,l)=dd; 

DEC(ke,2)=6; 

elseif ((aa==2)&(bb==2)&(cc==2)) 
DEC(ke,l)=dd; 

DEC(ke,2)=7; 

elseif ((aa== 1 )&(bb==2)&(cc==2)) 
DEC(ke,l)=dd; 

DEC(ke,2)=8; 

elseif ((aa== 1 )&(bb== 1 )&(cc==2)) 
DEC(ke,l)=dd; 

DEC(ke,2)=9; 

elseif ((aa== 1 )&(bb== 1 )&(cc==3)) 
DEC(ke,l)=dd; 

DEC(ke,2)=10; 

elseif ((aa==0)&(bb==l)&(cc==3)) 
DEC(kc,l )=dd; 

DEC(ke.2)=l 1; 

elseif ((aa==0)&(bb==0)&(cc==3)) 
DECike.l )=dd; 

DEC(ke.2)=12; 

elseif ((aa==0)&(bb==0)&(cc==4)) 
DEC(ke,l)=dd; 
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DEC(ke,2)=13; 

elseif ((aa== 1 )&(bb== 0)&(cc==4)) 
DEC(ke,l)=dd; 

DEC(ke,2)=14; 

elseif ((aa==l)&(bb==l)&(cc==4)) 
DEC(ke,l)=dd; 

DEC(ke,2)=15; 

elseif ((aa==l)&(bb==l)&(cc==5)) 
DEC(ke,l)=dd; 

DEC(ke,2)=16; 

elseif ((aa==2)&(bb==l)&(cc==5)) 
DEC(ke,l)=dd; 

DEC(ke,2)=17; 

elseif ((aa==2)&(bb==2)&(cc==5)) 
DEC(ke,l)=dd; 

DEC(ke,2)=18; 

elseif ((aa==2)&(bb==2)&(cc==4)) 
DEC(ke,l)=dd; 

DEC(ke,2)=19; 

elseif ((aa==3)&(bb==2)&(cc==4)) 
DEC(ke,l)=dd; 

DEC(ke,2)=20; 

elseif ((aa==3)&(bb==3)&(cc==4)) 
DEC(ke,l)=dd; 

DEC(ke,2)=2 1 ; 

elseif ((aa==3)&(bb==3)&(cc==3)) 
DEC(ke,l)=dd; 

DEC(ke,2)=22; 

elseif ((aa==2)&(bb==3)&(cc==3)) 



DEC(ke,l)=dd; 

DEC(ke,2)=23; 

elseif ((aa==2)&(bb==4)&(cc==3)) 
DEC(ke,l)=dd; 

DEC(ke,2)=24; 

elseif ((aa==2)&(bb==4)&(cc==2)) 
DEC(ke,l)=dd; 

DEC(ke,2)=25; 

elseif ((aa==l)&(bb==4)&(cc==2)) 
DEC(ke,l)=dd; 

DEC(ke,2)=26; 

elseif ((aa== 1 )&(bb==3)&(cc==2)) 
DEC(ke,l)=dd; 

DEC(ke,2)=27; 

elseif ((aa== 1 )&(bb==3 )&(cc== 1 )) 
DEC(ke,l)=dd; 

DEC(ke,2)=28; 

elseif ((aa==0)&(bb==3)&(cc==l)) 
DEC(ke,l)=dd; 

DEC(ke,2)=29; 

elseif ((aa==0)&(bb==2)&(cc== 1 )) 
DEC(ke,l)=dd; 

DEC(ke,2)=30; 

elseif ((aa==0)&(bb==2)&(cc==0)) 
DEC(ke,l)=dd; 

DEC(ke,2)=3 1 ; 

elseif ((aa==l )&(bb==2)&(cc==0)) 
DEC(ke,l)=dd; 

DEC(ke,2)=32; 
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elseif ((aa== 1 )&(bb== 1 )&(cc=0)) 
DEC(ke,l)=dd; 

DEC(ke,2)=33; 

elseif ((aa== 1 )&(bb= 1 )&(cc=l )) 
DEC(ke,l)=dd; 

DEC(ke,2)=34; 

elseif ((aa==2)&(bb== 1 )&(cc==l )) 
DEC(ke,l)=dd; 

DEC(ke,2)=35; 

elseif ((aa==2)&(bb=0)&(cc==l)) 
DEC(ke,l)=dd; 

DEC(ke,2)=36; 

elseif ((aa==2)&(bb==0)&(cc=2)) 
DEC(ke, 1 )=dd; 

DEC(ke,2)=37; 

elseif ((aa==3)&(bb==0)&(cc=2)) 
DEC(ke,l)=dd; 

DEC(ke,2)=38; 

elseif ((aa==3)&(bb=l)&(cc=2)) 
DEC(ke,l)=dd; 

DEC(ke,2)=39; 

elseif ((aa=3)&(bb== 1 )&(cc==3)) 
DEC(ke, 1 )=dd; 

DEC(ke,2)=40; 

elseif ((aa==2)&(bb==l)&(cc=3)) 
DEC(ke,l)=dd; 

DEC(ke,2)=4 1 ; 

elseif ((aa==2)&(bb==2)&(cc=3)) 
DEC(ke,l)=dd; 
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DEC(ke,2)=42; 

elseif ((aa==2)&(bb==2)&(cc==4)) 
DEC(ke,l)=dd; 

DEC(ke,2)=19; 

elseif ((aa==l)&(bb==2)&(cc==4)) 
DEC(ke,l)=dd; 

DEC(ke,2)=-15; 

elseif ((aa== 1 )&(bb==3)&(cc==4)) 
DEC(ke,l)=dd; 

DEC(ke,2)=-20; 

elseif ((aa== 1 )&(bb==3)&(cc==5)) 
DEC(ke,l)=dd; 

DEC(ke,2)=-19; 

elseif ((aa==0)&(bb==3)&(cc==5)) 
DEC(ke, 1 )=dd; 

DEC(ke,2)=-18; 

elseif ((aa==0)&(bb==4)&(cc==5)) 
DEC(ke,l)=dd; 

DEC(ke,2)=43; 

elseif ((aa==0)&(bb==4)&(cc==4)) 
DEC(ke,l)=dd; 

DEC(ke,2)=44; 

elseif ((aa==l )&(bb==4)&(cc==4)) 
DEC(ke,l)=dd; 

DEC(ke,2)=-21; 

elseif ((aa==l)&(bb==3)&(cc==4)) 
DEC(ke,l)=dd; 

DEC(ke,2)=-20; 

elseif ((aa==l )&(bb==3)&(cc==3)) 
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DEC(ke, 1 )=dd; 

DEC(ke,2)=-8; 

elseif ((aa==2)&(bb==3 )&(cc==3 )) 
DEC(ke,l)=dd; 

DEC(ke,2)=23; 

elseif ((aa==2)&(bb==2)&(cc=3)) 
DEC(ke,l)=dd; 

DEC(ke,2)=42; 

elseif ((aa==2)&(bb==2)&(cc==2)) 
DEC(ke,l)=dd; 

DEC(ke,2)=7; 

elseif ((aa==3)&(bb==2)&(cc==2)) 
DEC(ke,l)=dd; 

DEC(ke,2)=-24; 

elseif ((aa==3)&(bb==l )&(cc==2)) 
DEC(ke,l)=dd; 

DEC(ke,2)=39; 

elseif ((aa=3)&(bb== 1 )&(cc== 1 )) 
DEC(ke,l)=dd; 

DEC(ke,2)=-13; 

elseif ((aa=2)&(bb=l)&(cc=l)) 
DEC(ke,l)=dd; 

DEC(ke,2)=35; 

elseif ((aa==2)&(bb==0)&(cc==l)) 
DEC(ke,l)=dd; 

DEC(ke,2)=36; 

elseif ((aa==2)&(bb==0)&(cc==0)) 
DEC(ke, 1 )=dd; 

DEC(ke,2)=45; 
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elseif ((aa==l)&(bb==0)&(cc==0)) 
DEC(ke,l)=dd; 

DEC(ke,2)=-25; 

elseif ((aa== 1 )&(bb== 1 )&(cc==0)) 
DEC(ke,l)=dd; 

DEC(ke,2)=33; 

elseif ((aa==l)&(bb==l)&(cc==l)) 
DEC(ke,l)=dd; 

DEC(ke,2)=34; 

elseif ((aa==0)&(bb== 1 )&(cc== 1 )) 
DEC(ke,l)=dd; 

DEC(ke,2)=46; 

elseif ((aa==0)&(bb==2)&(cc== 1 )) 
DEC(ke,l)=dd; 

DEC(ke,2)=30; 

elseif ((aa==0)&(bb==2)&(cc==2)) 
DEC(ke,l)=dd; 

DEC(ke,2)=-2; 

elseif ((aa== 1 )&(bb==2)&(cc==2)) 
DEC(ke,l)=dd; 

DEC(ke,2)=8; 

elseif ((aa== 1 )&(bb==3)&(cc==2)) 
DEC(ke,l)=dd; 

DEC(ke,2)=27; 

elseif ((aa==l)&(bb==3)&(cc==3)) 
DEC(ke,l)=dd; 

DEC(ke,2)=-8; 

elseif ((aa==2)&(bb==3)&(cc==3)) 
DEC(ke,l)=dd; 
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DEC(ke,2)=23; 

end 

ke=ke+ 1 ; 
end 
end 

% Plot the transfer function of the 3-channel RSNS 
figure(l) 

plot(DEC( : , 1 ),DEC( : ,2)), grid 

xlabel('ADC INPUT V OLT AGE'),ylabel(' ADC DECIMAL OUTPUT CODE') 
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